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Abstract 

We study the radiative decays B — > V7 (where V = p,K*) in the framework of QCD fac- 
torisation, and in particular the hard-spectator contributions to the decay amplitudes. For 
the phenomenologically significant chromomagnetic operator, we show by an explicit next-to- 
leading-order computation that the spectator interactions factorise in the heavy-quark limit, i.e. 
that they can be written as the convolution of a hard-scattering kernel, computable in perturba- 
tion theory, and of the light-cone distribution amplitudes for the B and V mesons which contain 
the soft physics. The presence of an intermediate scale of O(MbAqcd) leads to the presence 
of Sudakov logarithms. We indicate how the demonstration of factorisation can be extended to 
other (four-quark) operators. 
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1 Introduction 



Radiative B — ► ^7 decays (where V is a light vector meson) are processes of particular interest 
in flavour physics which are already accessible at S-factories. Current measurements yield the 
following branching ratios [1]: 

B(B° -► K*°j) = (4.18 ± 0.23) • 10~ 5 B(B + -► iT + 7 ) = (4.14 ± 0.33) • 10~ 5 , 

whereas only upper bounds are available for the w°7, p°7 and /? + 7 modes (of order O(10 -6 ) at 
90 % confidence level) [2]. Within the standard model these decays can give measurements of 
the V-td and V* s elements of the CKM matrix. 

Penguin- mediated processes, such as b — > 57 and b — > radiative decays, exhibit particular 
sensitivity to physics beyond the standard model. It is therefore theoretically important to 
be able to disentangle the Standard Model contributions from potential effects of new physics. 
The principal difficulty in calculating the amplitudes is due to the presence of non-pertubative 
QCD effects. For inclusive decays, such as B — > X s ^f, an operator product expansion allows the 
decay amplitude to be computed as a series in inverse powers of TO5 (the mass of the 6-quark), 
with the leading term being calculable in perturbation theory [3]. For exclusive decays the 
decay amplitudes can be calculated in the framework of QCD factorisation, first developed for 
nonleptonic two-body decays [4, 5]. At leading order in 1/Mb, the long-distance QCD effects are 
factorised into universal quantities, the B — > meson transition form factors and the light-cone 
distribution amplitudes of the mesons. 

The factorisation framework has been applied to B — > K*^y and B — > P7 decays, refs. [6, 7], where 
in addition to the leading-twist contributions, the subleading (but phenomenologically relevant) 
effects of weak annihilation were also considered. The same class of decays was also analysed 
in an effective-theory framework at higher orders in perturbation theory [8]. The subsequent 
phenomenological analysis suggests interesting quantities to constrain the shape of the unitarity 
triangle, such as the CP-asymmetry for the 7 angle [6] and the ratio of branching ratios 
for p°7 and K*°~f for the Rt side (which provides bounds already competitive with the B — B 
and B s — B s mass differences) [9] . These flavour-changing neutral-current processes yield strong 
constraints on supersymmetric models [10]. 

In view of their phenomenological importance, B — ► V~i decays deserve further investigation in 
the framework of QCD factorisation. In this paper we focus on (strong) radiative corrections 
to spectator interactions, in order to determine whether and how factorisation holds at higher 
orders in perturbation theory. We have previously followed a similar line of investigation for the 
purely radiative decay B — > ^tv [11] (and the related decays B — > 77 and B — > ^£ + £~ [12]) in 
order to acquire a better understanding of the properties of the light-cone distribution amplitude 
of the .B-meson. B — > V7 decays provide an opportunity to study similar issues concerning the 
factorisation of long-distance effects, but now with the presence of a hadron in the final state. 
In addition, these decays possess similar features to heavy-to-light semileptonic decays, such as 
B — > tt£u which are currently being studied in the SCET framework [13, 14] (see refs. [15, 16, 17] 
for recent progress). 

The effective Hamiltonian for B — > Vj decays, where V is a light-vector meson V is [18]: 




(1) 
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where X^ = V* q V p b (unitarity of the CKM matrix implies that = — + A^) and so 
contributions from diagrams with loops containing a top quark are included implicitly). Q\ - 
Qq are four-quark operators: 

QP = (qp) v ^ A (pb) V -A Q\ = {q~iPj)v~A (pj h)v-A (2) 

Qs = (qb)v-A ^(q'q')v-A Qa = (ft bj) V -A ^(<?j q'i)v-A (3) 

q> q< 

Q 5 = (qb) V -A ^{q'q')v+A Q% = (ft bj) V -A ^{q- q'i)v+A (4) 

q> q< 

and Q7 and Qs are the electromagnetic and chromomagnetic penguin operators: 

& = -^rn b qa^(l + 75)6 ^ Qs = ^m b qa^(l + 75 )T A b G% . (5) 

q = d or s and the convention for the sign of the coupling constants corresponds to the covariant 
derivative = 8^ + ieQfA^ + igT a A^, with and A^ representing the photon and gluon 
fields respectively and Q e = — 1 etc. The Wilson coefficients Cj are known at next-to- leading 
order in renormalization group improved perturbation theory [19]. 

The calculation of the amplitude for B — > V^7 decays requires the evaluation of the matrix 
elements of the weak operators Qi listed above. We will restrict our analysis to the phenomeno- 
logically relevant operators Qi, Q7, Qs'- the QCD penguins Q3 . . . Qe start at 0{a s ) and are 
multiplied by small Wilson coefficients in the weak Hamiltonian TC, whereas the contribution 
from Q2 starts only at 0(a 2 s ). 

The factorisation formula for B — > Vj decays is of the form 



(6) 



where the non-perturbative effects are contained in F B ^ v ,a form factor for B — > V transitions, 
and $ B and the leading- twist light-cone distribution amplitudes of the B and V mesons. 
The hard-scattering amplitudes T/ and T/ 7 include only short-distance effects and are calculable 
in perturbation theory. The indices p and q have been suppressed on the right-hand side of 
eq. (6). It is expected that the factorisation formula is valid up to corrections of 0(AQCD/?rift)- 
Following the analysis in refs. [6, 7], the contribution from the electromagnetic operator Qj is 
included in the first term on the right-hand side of eq. (6), corresponding to a form factor. The 
chromomagnetic operator Qs and the four-quark operator Q\ contribute to both terms. 

The subject of this paper is the second term on the right-hand side of eq. (6), which arises 
from hard spectator interactions. For such terms factorisation has been explicitly established 
only at leading order in perturbation theory. We consider the spectator contribution to the 
matrix element of the chromomagnetic operator Qs, and demonstrate by explicit calculation 
that the mass singularities at one-loop order are precisely those of the distribution amplitudes 
in eq. (6). The detailed results for the mass singularities from each diagram are presented in 
appendix A, together with the terms containing "large logarithms". We analyse the reasons 
for the factorisation of mass singularities, and present a heuristic demonstration for the matrix 
elements of Qs in section 3. We then argue that such a cancellation is also valid for the matrix 
elements of the remaining operators (appendix B), analysing the contribution from Q\, which is 
expected to be the largest from the four-quark operators Q\ - Qq, in particular detail in section 4. 
Section 5 contains our conclusions. 



2 



1.1 Kinematics 



We consider the decay 

B(p)^ 7 ( qi ,e*)V(q 2 , V *) (7) 

where p, q± and q 2 are the momenta of the corresponding particles and e* and rj* are the 
polarization vectors of the photon and V respectively. We work in the rest frame of the B- 
meson, p = (Mb,0), and take the momenta of the photon and V to be in the z-direction, with 
both final-state vectors transversely polarized. 

It is convenient to introduce light-cone coordinates; for any four- vector I we write I = (l + ,l-,l±), 
with 

l± = go that l 2 = 21+1- - ll . (8) 

We neglect the mass of the vector meson, and take q\ and q 2 to be in the + and — directions 
respectively: 

« = M' ?2= (°'7¥' 51 )- <9) 



1.2 Evaluating the T/ 7 

In this paper we focus on the contribution from the hard spectator interactions represented 
by the second term on the right-hand side of eq. (6), which we write in the schematic form 
<fr B ® Tj ® Here and in the remainder of the paper, for simplicity of notation, we suppress 
the superscript II. 

In order to verify the validity of the factorisation formula (6) we can use perturbation theory 
with conveniently chosen partonic external states. In particular we evaluate the hard-scattering 
kernels Tj , and a signal of the breakdown of factorisation would be the presence of residual long- 
distance effects in the Tj. For our calculations we take the initial state to have momentum p and 
to consist of a 6-quark with momentum p — k and a light (spectator) antiquark q' with momentum 
k. The components of k are of O(Aqcd)- Our final state consists of the photon with momentum 
qi and a quark (q) and spectator antiquark (q 1 ) with momenta xq 2 and xq 2 respectively. The 
partons in the initial and final states are on their mass shells. The hard-scattering kernel is 
obtained in the standard way; by using perturbation theory to calculate the matrix elements 
of the operator Qj between these states together with the light-cone distribution amplitudes 
(defined explicitly in the next subsection) for the initial and final partonic states. Let Ai denote 
the left hand side of eq. (6) and the superscript (n) label the perturbative contributions of 0{a™) 
relative to the lowest order (denoted by (0)). The hard-scattering kernels at one-loop order (T^) 
are then obtained by rewriting the factorisation formula as 

(o) T (i) (o) = ^(l) _ $6 <r (i) T (o) (o) _ (o) r (o) (i) (10) 

where the $'s denote the corresponding distribution amplitudes. We confirm below that al- 
though there are mass singularities in the they are cancelled by the remaining terms on 
the right-hand side of eq. (10). 

Throughout this paper the discussion of the contributions from particular diagrams corresponds 
to the Feynman gauge. 
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1.3 Light-Cone Distribution Amplitudes 

We define the light-cone distribution amplitude of a state H containing the 6-quark by 

= / dz ~ e fk + z -(0\^(z)[z,0}b a (0)\H)\ z+ , z±=0 , (11) 

where q',b are the quark fields and a, (3 are spinor labels (we denote the spectator antiquark by 
q'). [z,Q] denotes the path-ordered exponential, "Pexpf— ig j Q dx fJ 'A^(x)]. Ultimately of course, 
when evaluating the physical decay amplitude we use the distribution amplitude of the S-meson, 
but in the evaluation of the hard-scattering amplitude T; L in perturbation theory we take H to 
be a quark-antiquark state. 

Similarly the light-cone distribution amplitude of a light state L is defined by: 

*£(u) = ft _| % ~/ )+ e -^+^) (Ll^^^xl^^lO)!^).,^)^ , (12) 
where u = 1 — u, and the integrand in eq. (12) is a function of (x — y) + . 

In terms of the definitions in equations (11) and (12), for the B and V mesons we follow an 
equivalent spinor decomposition to that in refs. [20, 5, 6] and define the distribution amplitudes 
$±(fc+) and by 2 

= -^^{(l + ^[^(fc + ) + ^(^(^ + )-$^(fc + ))] 75 }^ (13) 

where the light-like vector n_ is given in light-cone coordinates by n_ = (0, \/2, OjJ, and the 
four velocity v is defined by = Mbv^. The decay constants fs and fy are defined by 

(0\q' W b\B(p))=if BPfl and r?*) j^g' | 0) = -i^K - tfM ft . (15) 



2 Leading Order Results for the Matrix Elements of the Oper- 
ators Q 8 and Qi 

Before studying the cancellation of mass singularities at next-to-leading order (NLO) in pertur- 
bation theory, we present the lowest-order contributions to the hard-scattering amplitude from 
the operators Qs and Qi respectively. 



2.1 Tree-Level Result for the matrix element of Q 8 



As explained above, we determine the hard-scattering amplitude Tg by choosing the convenient 
external states H = bq' and L = qq' and computing both the matrix element 

As^-^^iL^QslH) (16) 

2 Following ref. [11, 12], $± is defined in the Heavy-Quark Effective Theory, but with the physical decay 
constant fs in front of it. It was pointed out in ref. [21] that a different normalisation factor can be chosen to 
ensure a more natural comparison with HQET results. Such a change in the normalisation of does not affect 
the outcome of our analysis and corresponds to a simple redistribution of terms between the hard-scattering kernel 
T 11 and the B-meson's light-cone distribution amplitude. 
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Figure 1: Lowest-order contributions to the matrix element (q'{xq2) q(xq 2 )"f(qi) I Qs I l'(k) b(p- 
k)). The spring-like line represents the photon and the square denotes the insertion of Q$- 



and the distribution amplitudes <& H ' L . In this section the calculation is performed at lowest 
order of perturbation theory. 

There are four diagrams contributing to at lowest order, and these are shown in fig. 1. At 
leading twist only the diagram of fig. 1(a) contributes and we find: 

4 0) = -G F \^ ^ _L (u(xq2){flu(l + l 5 )T A }u( P -k)) (v(k) {~fT A } v{xq 2 )) , (17) 

where u and v represent the free-particle spinor wave functions and Q q = —1/3 is the charge of 
the quark q. 

The distribution amplitudes at lowest order are: 

$ S (0) (M = 2ir5(k + - ~k + ) vp(k)u a (p - k) (18) 
*$ (0) («) = - x) u s (xq 2 ) Vl {x q2 ) , (19) 

and can be written in the factorised form 

^(o) = $ 6 9 -'(o)^ r (%$^(0) (20) 

with 

T^s(h,u) = -G F X? ^ 4- {A,(l + 7s)T A } 5a {l v T A } M . (21) 

For simplicity of notation, we have suppressed the colour indices. is a function of the 

convolution variables k + and u, but does not depend on any kinematical quantities of 0(Aqce>)- 
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Figure 2: Lowest-order contribution to the matrix element {q'{xq2) q{xq2) \ Qi \ q'{k) b{p — k)). 
The arrow on the quark loop defines the direction of the flow of fermion number. 



We conclude the lowest-order calculation by briefly considering the three diagrams in fig. 1 (b) - 
(d). Fig. 1(b) manifestly gives a higher- twist contribution, since both the gluon and internal 
antiquark are off-shell by an amount of O(Mg). In diagram 1(c), although the gluon is off- 
shell by O(Mg), the internal antiquark is off-shell by 0((q\ — k) 2 ) = 0(M#Aqcd)- However 
this potentially leading-twist contribution vanishes when projected with the wave function of a 
transversely polarized vector meson (see eq. (14) ). The contribution from diagram 1(d) vanishes 
by the equations of motion of the final-state meson. 



2.2 Lowest-order result for the matrix element of Qi 

The second example which we study in some detail is the matrix element of the operator Q\. 
At lowest order the relevant diagrams are shown in fig. 2. The contribution proportional to 
(Xu ^) is the difference of the diagrams with the charm (up) quark propagating in the loop and 
one with the top quark. Here we present the results for the charm quark (the corresponding 
results for the up quark can readily be deduced from these): 

Al = mV ~ A {db)v - A lH) ' (22) 

Following ref . [6] , we do not neglect m 2 / m 2 , and hence keep an important imaginary part which 
vanishes when m c is set to zero. 



Evaluating the diagrams in fig. 2 we find: 

4 0) = eQ £ S xq2 \ + (u(xqMl- l5 )T a u(p-k)) (v(kh P T a v(xq 2 

{T 1 (x,s)e^ Xu (xq2-qi)a + T 2 (x,s)e^ TV q p 1 q la q 2T + T 3 (x,s)e aprX q u 1 qi^2r} X e* x , (23) 

where Q c = 2/3 is the charge of the charm quark, s = m 2 /m 2 , 

1 4s 2s 

T^s) = -- + L, (24) 

3 x x 

m , x 1 /2 8s 4s \ 



4 8s 4s 



T ^- s) - -w B {s + --- (L + p >i' < 26 > 



with 



L = V 1 — 4s/x + le log 

[^l-As/x + ie-l 



(27) 



6 



and 

7== J + ^ f 7^—1 - (28) 

where L 2 is the dilogarithm function. From eq. (23) and the tree-level expressions for the 
distribution amplitudes in equations (18) and (19) one immediately obtains the lowest-order 
contribution to the hard-scattering kernel T^°\ 

2.3 Predictions for the Spectator Contributions to B — > Vj Decay Ampli- 
tudes at Lowest Order 

In order to obtain the predictions for the spectator contributions to B — > V^f decay amplitudes, 
the hard-scattering kernels obtained in the previous two subsections need to be convoluted with 
the B and V distribution amplitudes. Convoluting the lowest-order result for the hard-scattering 
kernel Tg "* in eq. (21) with the distribution amplitudes we obtain for the spectator contribution 
to the amplitude 

.f^Wifi) - - (Gf a««,^m^ {eQq) jtkim j du m 

x {e^e^rj^q^ + iie* ■r j *)(q 1 -q 2 )] , (29) 

where Q q = —1/3 is the charge of the down or strange quark and N = 3 is the number of 
colours. This is equivalent to the result in eqs. (39) and (40) of ref. [6]. 

Similarly, convoluting the hard-scattering kernel from the lowest-order amplitude in eq. (23) 
with the B and V-meson distribution amplitudes we obtain 



g fW n, \nc\R\ G^a s C F / B /^ f dk + ^fc+) f , 

-^(VimB) = — — — (eQ^j—^—jduZyMHu,*) 

x V%V + ■ V*)(Qi ■ <fc)} , (30) 



where 



h(u, s) = -- + ^ { L 2 I p I + L 2 



II II- I \ 1 _ ^j u-is+ie J \ 1 + ^/ u ~ 4s+ie 



(31) 



If m c is set to zero, h(u, 0) = —2/u. These results are equivalent to those in eqs. (33) -(35) of 
ref. [6]. 



3 Factorisation at One-Loop Order for the Chromomagnetic 
Operator Q 8 

We now turn to the main subject of this paper, the factorisation of mass singularities at one- 
loop order for the spectator interactions. We start with the operator Qs and study the one-loop 
corrections to the lowest-order diagram in fig. 1(a). We redraw this diagram in fig. 3, labelling 
the external lines by (1)- (4), the internal propagators by (5) and (6) and the weak vertex by 
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Figure 3: Notation for one-loop diagrams contributing to the B — * Vj decay from Qg. The 
extra gluon can be attached to any pair of circles (which might be identical). 



(7). The mass singularities at NLO arise from diagrams in which a gluon is attached to two 
of the circles in fig. 3 (including diagrams in which both ends of the gluon are attached to the 
same circle). We denote by (ij) the diagram obtained by adding a gluon between the circles (i) 
and (j) (see, for example, fig. 4 where diagrams (34), (23) and (13) are drawn explicitly). We 
evaluate the terms containing the mass singularities and/or large logarithms which are of the 
form log(Mg/<72.A;). 

We present the results from an explicit evaluation of each diagram in Appendix A, from which the 
cancellation of mass singularities from the hard-scattering kernel is apparent. In the calculation 
presented in the appendix we do not make any a priori assumptions about the regions of phase 
space which lead to mass singularities or large logarithms. In this section we demonstrate the 
cancellation of mass singularities using a heuristic argument, based on power counting, which 
shows that the mass singularities arise from regions in which the gluon is soft or collinear 
with q2 or k (including contributions from the soft- collinear region of phase-space [24]). In 
Appendix B we argue further that such a cancellation holds for all the weak operators in the 
effective Hamiltonian of eq. (1). 

We emphasize that not all the one- loop corrections are given by the set of diagrams {(ij)}, and 
for illustration in figure 9 we draw an additional Feynman diagram (called D). However, all 
the mass singularities are given by diagrams (ij) and only these diagrams and diagram D give 
large logarithms. The remaining diagrams, which may contribute either to T l or T n , are only 
sensitive to physics at scales of 0(Mb)- 

We implement a similar notation for diagrams which give the one-loop corrections to the dis- 
tribution amplitudes, and which therefore contribute to <3? b5 ' W (g> ® <& 9 ^ (o) an( j $6?' (o) ^ 
$ $99' W (see eq. (10) ). By (13)', (11)' and (33)' we denote the contribution to (i) g, 
J 1 ' ' (g> <3?' ? ^ (°) in which the gluon in the one-loop correction to the distribution amplitude of the 
initial bq' state (cJ^'W) j s attached to lines 1 and 3 as indicated (where 1 denotes the 6-quark 
and 3 the spectator antiquark). Similarly by (24)', (22)' and (44)' we denote the contribution to 
$bq' (0) ^ T (o) <g, $<?<?' (i) j n w hich the gluon in the one-loop correction to the distribution amplitude 
of the final qq' state (3> 99 " is attached to lines 2 and 4 as indicated (where 2 and 4 denote 
the quark and spectator antiquark respectively). There are also corrections to the distribution 
amplitudes in which one end of the gluon is attached to the path-ordered exponential. By (IB)' 
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{(SB)') we denote the contribution to U) $ ^(0) ^ (0) in w hi c h one end of the gluon in 
the one-loop correction to the distribution amplitude of the initial bq' state is attached to line 1 
(line 3) and the other to the path-ordered exponential. Similarly by (2V)' ((4V)') we denote the 
contribution to <°) (g> T< ) <g> W in which one end of the gluon in the one-loop correction 
to the distribution amplitude of the final qq' state is attached to line 2 (line 4) and the other to 
the path-ordered exponential. 

We now present our heuristic diagrammatic argument for the cancellation of the mass singu- 
larities. In Appendix B we give a formal argument explaining the reasons for this cancellation 
based on the collinear Ward identity. 



3.1 Cancellation of Soft Divergences 

We start by considering soft divergences, which arise when the loop momentum is small, with all 
components vanishing in a similar way. Such divergences arise in diagrams with a gluon attached 
to two external lines. As explained below, some of these divergences are absorbed into the 
mesons' distribution amplitudes in a straightforward manner. The cancellation of the remaining 
infrared divergences occurs as expected from arguments based on colour transparency [22] . Soft 
gluons only couple to the compact, colour-singlet vector state through a dipole interaction, which 
is non-singular. In terms of Feynman diagrams there is a cancellation of infrared divergences 
between graphs in which the soft gluons couple to different constituents of the vector meson. 

To illustrate this point consider the graph (34), drawn explicitly in fig. 4, which is singular in 
the soft region, V 3, — > 0. The contribution from this region is 

*(n C A a(P) f dH ±(k-l)-xq 2 



(34) = -ig z \C F ^ A 



2 ) ~ 8 J (2ir)*l 2 (l-k)i(xq 2 -l) 2 
~ 2ia 2 (cr-^S A {0) f dH {k ~ l) ' Q2 (32) 

In the soft region the contribution from diagram (23) (also drawn in fig. 4) is 

i( n Ca \ a(°) f 4:{k-l)-xq 2 



(23) = ig*lC F --£)A 



2 J 8 J (2ir)U 2 (l-k) 2 (xq 2 -l) 2 
-»f (O - ^) 4 0) / T~r~TT r','/' H:!'' 1 „ = "(34). (33) 



' f >• _ 9A\ a(°) [ ^ ( fc ~ • g2 

2 J^ 8 J (27r) 4 P(l-k) 2 (q 2 .l) 

The soft divergences in diagrams (12) and (14) cancel in an analogous way. In addition there 
are soft divergences in the self-energy diagrams (nn) (n=l-4), but these are trivially absorbed 
into the B and vector-meson light-cone wave functions (diagrams (nn) 1 , with n=l-4). 

Before proceeding to a discussion of collinear divergences, we briefly comment on diagram 
(13) drawn in fig. 4. By power counting one can readily demonstrate that this diagram has 
a leading- twist contribution from the soft region in which the components of I are O(Aqcd)- 
This contribution is also absorbed into the corresponding distribution amplitude; specifically 
the contribution from diagram (13) is cancelled by the subtraction of the corresponding term in 
$bg' (i) ^ t(o) ^ ^qq' (o) ( Le _ by the SUD traction of dia gram (13)' in which the diagram contribut- 
ing to W is also drawn in fig. 4). The reason for such cancellations has been explained in 
detail in ref. [11]. Diagram (13)' differs from (13) in that the propagator of the gluon attached 
to the weak vertex is replaced by its eikonal approximation ((xq 2 — k + l) 2 — > —2xq 2 ■ (k — I)). 
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z 



(13) (13)' 

Figure 4: The one-loop graphs (34), (23) and (13) contributing to the decay amplitude and a 
diagram which contributes to the distribution amplitude of the initial state W . The dashed 
line represents the path-ordered exponential and z 2 = 0. 

Since the leading-twist contribution in (13) comes from the soft region, we can approximate 
this propagator by the eikonal form in this diagram. Thus the two contributions are equal and 
therefore cancel in the evaluation of the hard-scattering amplitude. Diagrams (24) and (24)' do 
not have such leading-twist contributions from the soft region. 

3.2 Cancellation of Collinear Divergences 

We now consider collinear divergences, present when a gluon is attached to an external light- 
quark line. Such singularities can occur when the gluon is parallel to the momentum of the 
vector meson q 2 or to the momentum of the spectator antiquark k (in this latter case, since 
k = O(Aqcd); one does not call the corresponding gluons collinear in the standard SCET 
nomenclature). We distinguish between these two cases, beginning our discussion with the 
former. 

3.2.1 Divergences from the Region Collinear to q-z 

In this section we consider the region of phase space in which the loop momentum is parallel 
to q 2 . There are corresponding mass singularities in diagrams with a gluon attached to one of 
the final-state external lines, 2 or 4. Before evaluating each diagram in turn, it is convenient to 
consider the distribution amplitude of the final state, in this case the state with a light quark 
and a light anti-quark with momenta xq 2 and xq2 respectively: 

®%{u)=q 2 -f ^±e-^ x+ ^ (q'(xq 2 )q(xq 2 )\q^y)[y,x}q a (x)\0) , (34) 
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y xq 2 y xq 2 p -k o 

EZ 

(a) (b) (c) 

Figure 5: Two one- loop diagrams contributing to the distribution amplitude of the qq' final 
state, (a) and (b), and a diagram contributing to the distribution amplitude of the bq' final 
state. The dashed line represents the path-ordered exponential and z 2 = (x — y) 2 = 0. 



where z = x — y, a and j3 are spinor labels and [y, x] represents the path-ordered exponential 
between the two points y and x. 

We have seen above that for the free theory (see eq. (19) ) 

®% \ u ) = 5{u-x) u/3{xq 2 ) v a {xq 2 ) , (35) 

where ur and v a are free Dirac spinors. Now consider the contributions to the distribution 
amplitude from the two one-loop diagrams in fig. 5(a) and (b), in which a gluon is attached at 
one end to the path-ordered exponential 3 . In the collinear region in which I, the momentum of 
the gluon, is parallel to q 2 (I — £q 2 , with the remaining components of I vanishing as Z+ ~ A 2 /Z_ 
and l± ~ A as A — > 0) their contributions can be written as: 

^SSW = ^(^2) v a (xq 2 ) x 

(V C F) J i^-L-p {6{u _,)_,(„_,_ {)} , (36) 



= up{xq 2 )v a {xq 2 ) x 

° f) I T*T^hw {6{U - X ^-^ U - X + , (37) 

where the subscripts (a) and (b) refer to the corresponding diagrams in fig. 5. We will use the 
representations in eqs. (36) and (37) in our discussion below. For the other diagrams contributing 
to the wave function of the final state at one-loop order it will be very straightforward to see 
how the mass singularities are precisely those needed to absorb the corresponding ones from the 
amplitude. 

In order to illustrate the cancellation of collinear divergences coming from the region of phase 
space in which the momentum of gluon(s) in the loop is(are) parallel to q 2 consider diagram (12). 
Since we are considering the collinear divergence we set / ~ £q 2 , with £ finite and the remaining 
components of I vanish as /+ ~ A 2 //_ and l± ~ A in the singular region A — > 0. Keeping only 
the leading-twist terms one readily finds 

The integral over /+ and l± gives a collinear divergence. These collinear divergences also cancel, 
but in a different way to the soft divergences. In this case it is straightforward to verify that in 
the region where / is collinear to q 2 , 

(23) = -(12), (39) 

3 The remaining one-loop diagrams contributing to the distribution amplitude will be discussed later. 
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so that the corresponding collinear divergences cancel between diagrams (12) and (23). As we 
will see when we consider the cancellation of collinear divergences for the operator Q\ below, 
such a cancellation between pairs of diagrams is not typical. What is required however, is that 
apart from the propagators which are explicitly exhibited in eq. (38) (l/l 2 and l/(x<72 — I) 2 ) one 
can replace Z by ^ everywhere. 

There are similar cancellations between other pairs of diagrams 4 

<">=-<«> - -^-fMi^^CT (40) 

(26) = -(27) = ^4°>?/|LI_L_ ; (41) 
(36) = -(37) = -iiC A Af f ^j^. (42) 

This leaves us with the collinear divergences from diagrams (25) and (45) in which one end of 
the additional gluon is attached to the internal quark propagator 5. These are precisely the 
terms which are absorbed into the wave function of the final state (in this case the qq' state), 
i.e. the divergences in the contributions from diagrams (25) and (45) to are cancelled by 
the subtraction of <S> b «' (°) (g> T 8 (0) ® W in eq. (10). More specifically they are cancelled by the 
contributions to the wave function from the diagrams in fig. 5(a) and (b), which we denote by 
(2V)' and (4F)' respectively. Using the expressions in eqs. (36) and (37), it is straightforward 
to demonstrate that in this collinear region 

(25) = eyy—VftoAPf^j^; (43) 
„2^_ 4 (o) f dH g-j 1 1 



(45) = (wy = vf Cr JiS> jf^TTJp TT^TTF ' < 44 » 

Finally we discuss the contribution of diagram (24). For the processes considered in this paper, 
with a transversely polarized vector meson in the final state, the leading-twist contribution 
vanishes. We therefore do not discuss this further, other than to note that for a generic final 
state this contribution is cancelled by the corresponding one in 

Thus we find that, at one-loop order, there are no singular contributions to from the 

region in which the momentum of the gluon I is parallel to q2, and that the cancellation of the 
singularities has a very simple structure. The divergences either cancel between pairs of diagrams 
contributing to -4^, or between diagrams contributing to A^ and to <$> w (°) <g> T 8 (0) <g> W. 



3.2.2 Divergences from the Region Collinear to k 

For the 5-meson initial state, the momentum of the spectator antiquark, k, is generically soft, 
with all components of O(Aqcd)- In ° ur calculations we take k to be light-like so that long- 
distance effects are manifested by the presence of mass singularities, specifically collinear diver- 
gences from the region of phase space where the momenta of internal partons are collinear with 

4 The mass singularities considered in this section appear in diagrams in which a gluon in a loop is attached 
to one of the final-state external lines 2 or 4. Below however, we give singular expressions for diagrams (36) and 
(37). These diagrams could equally well have been denoted by (46) and (47) and hence are indeed singular. 
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k (as well as the soft divergences discussed above). We now consider these collinear divergences. 
The diagrams which can have such a divergence contain a gluon line attached to the external line 
3. Writing I ~ ok with Z_ and l± vanishing with A as A 2 // + and A respectively, we readily find 
the following three leading- twist contributions (diagram (37) =0 and (35) gives a higher- twist 
contribution) : 

(23 > . _ 2i/ ( CF _^) 4 o>/^I^_L_, (45) 

2 ( n C A f d 4 l 1-a 1 



(34) = V| C ,-^^y (46) 



(36) = -tfc^jgp^-y. (47) 
Summing these contributions we obtain 

(23) + (34) + (36) = -2 ig "C F J ^ p ^ = (35)', (48) 

where (3-B)' is the corresponding one-loop contribution to cJ^'W (gi (g $^'(°), i.e. the one 
from the contribution of the graph in fig. 5(c) to <J> b ^ Thus there are no collinear singularities 
in Tg 1 ^ arising from the region in which / is parallel to k. 



3.3 Double Logarithms 

In sections 3.1 and 3.2 we have analyzed the regions of phase space in which the momentum of 
a gluon in the loop was either soft or collinear with q2 or k, and deduced that the cancellation 
of the corresponding divergences is straightforward. We now briefly analyze the region in which 
the gluons are both soft and collinear, which leads to double logarithms in an infrared cut-off 
for some of the diagrams. 

The double logarithms arise in diagrams (12) and (14) as ^ — > (see equations (38) and (40)) 
and in diagrams (23) and (34) as either £ — ► (equations (38), (39) and (40) ) or a — > (equa- 
tions (45) and (46)). Explicit evaluation of the diagrams, described in Appendix A, confirms 
that only these four diagrams contain double logarithms. 

From the expressions in the previous subsections we can also deduce how the double logarithms 
cancel. We know from section 3.1 that the soft divergences cancel in the sum of diagrams 
(23) +(34), and from equations (45) and (46) we note that the term which diverges as 1/a as 
a — ► cancels in the integrand of (23) +(34). Thus the double logarithms in diagrams (23) and 
(34) cancel, leaving only a collinear divergence from the region in which / is parallel to q2- This 
is confirmed by explicit calculations (see eqs. (79) and (80) in Appendix A). Similarly, since the 
soft divergences cancel in diagrams (12) and (14), in the sum of these two diagrams we are only 
left with a single logarithm corresponding to a collinear divergence from the region in which I 
is parallel to c/2- 

When presenting such heuristic arguments, one has to consider the possibility that, after the 
cancellation of the double logarithms, a subtlety in the way that the integrals are regulated may 
leave a single logarithm which had not been considered. The explicit calculations in appendix 
A clearly remove such worries, but it is also relatively straightforward to verify this without the 
full computations. For example one may evaluate either the /_ or l + integration by contours, 
and readily verify the above pattern of cancellations. 
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p-k 




k 



® 



® 



Figure 6: Notation for the next-to-leading order diagrams contributing to the B — > Vj decay 
from Qi (the diagram with the photon emitted from the other propagator in the quark loop is 
included implicitly). The extra gluon can be attached to any pair of circles (which might be 
identical). (n7), for n = 1-6, represents the four diagrams for which the external photon and 
internal gluon are attached to either the quark or antiquark in the loop. 



4 Factorisation at One-Loop Order for the Four-Quark Opera- 
tor Qi 



In this section we repeat the heuristic arguments with Q\ as the weak operator. For the soft 
divergences which are present in diagrams in which the additional gluon is attached to two 
external lines, the cancellation occurs in the standard way, exactly as described in section 3.1 
and we don't discuss them any further. For the collinear divergences, we demonstrate below 
that they also cancel, but the pattern of cancellation is different in detail from that in section 3. 
Moreover the cancellation occurs even if we keep higher- twist contributions such as m^/m^ or 
2^2 • k/m^ in the quark- loop. This is consistent with the the arguments based on collinear Ward 
identities presented in appendix B. 

The notation for the next-to-leading order diagrams is similar to that in our discussion of the 
contributions from Qg. We denote by (ij) the diagrams in which the external photon is emitted 
from either of the two propagators in the quark loop and the internal gluon is attached to circles 
i and j. By (n7) we mean the sum of the diagrams with one end of the gluon attached to either 
the quark or antiquark in the loop. 

4.1 Divergences from the Region Collinear to q 2 
4.1.1 Neglecting terms of 0(m^/m^) 

We now repeat the steps of section 3.2.1 for the operator Q\ with the photon emitted from the 
internal loop, beginning with the simpler case in which we neglect m c . As before we start with 
diagram (12) (see fig. 7) for which the expression in this collinear region is now 



(12) = 2* 5 2 (C F - ^) 4 0) 



/ 



dH x-£ 1 1 



(49) 



(2vr)4 £ P (xq 2 - 1) 



2 ' 
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where A^ is the amplitude at tree level. This is no longer cancelled by the corresponding 
contribution from diagram (23) (see fig. 7) which is now given by 

Using eq. (36), we can see that the contribution (2V)' to (°) ® T g (0) ® cancels the 

abelian component of these two diagrams: 

<2^Vc F <>/!L^_L_ (51) 

Indeed the abelian component of diagram (12) is cancelled by the contribution to {2V) 1 coming 
from the first delta function in eq. (36) and similarly the collinear divergence in diagram (23) 
is cancelled by the contribution from the second delta function. This simple observation will be 
important when we generalise the discussion to include corrections of 0{m 2 /m 2 ). Thus we have 



(12) + (23) - (2V Y = -ii C A AT j^-^p • «*> 

In the case where we neglect the mass of the charm quark, the non-abelian contribution in 
eq. (52) is cancelled by diagram (26). There are no collinear divergences from the sum of 
diagrams (27). 

Similarly, the expressions for diagrams (34) (which is unchanged from that in eq. (40)), (14) and 
(41/)' (see eq. (37) ) are now 

(34 ) . W (c F -f)<>/|^l_L_ m 



d 4 l x 1 1 



(14) = -Vf^Cr-fJXPJwjtj^ry -d (54) 

W = -V Cf ^/|L^_L_, (55) 



so that the abelian component cancels and 

dH 



(U) + (34)-(4V)' = i 9 2 C A A? 



I 2 (xq 2 - If 



(56) 



This contribution is cancelled by diagram (46). If we neglect the mass of the charm quark there 
are no collinear divergences from the sum of the diagrams (47). 

Thus the collinear divergences from the region of phase space in which I is parallel to q 2 either 
cancel or are absorbed into the vector meson's distribution amplitude. 



4.1.2 Including terms of 0(m 2 ,/mf ) ) 

In this subsection we demonstrate that the collinear divergences cancel, even if we include the 
corrections of 0{m 2 c /m 2 ). The additional complication is that the integral over momentum in 
the charm-quark loop now depends on whether the outgoing momentum in the gluon attached 
to the loop is xq 2 — k (as in diagram (12) in fig. 7 for example) or xq 2 — k + I (as in diagram 
(23) in fig. 7 for example). We will consider diagrams (27) and (47), in which two gluons are 
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xq 2 



XQ2 



(12) (23) 

Figure 7: One-loop diagrams (12) and (23) contributing to corrections to the matrix element of 
Q\. The corresponding diagrams in which the photon is emitted from the other propagator in 
the quark loop are to be considered implicitly included. 




xq2 - k xq2 - k 

(a) (b) 
Figure 8: Two of the attachments of the gluon of momentum I onto the quark loop (see text) . 



attached to the charm-quark loop, separately. In order to distinguish between the two types of 
contribution we write 

Af\x) = -h(x), (57) 

x 

where the x dependence in I\{x) comes from the integral over the loop-momentum in the quark 
loop and is exhibited explicitly in eqs. (23) -(26). If we neglect m^/m^, then Ii(x) is indepen- 
dent of x, which leads to the vanishing of some contributions, such as (27) and (47). Since 
now we do not neglect these effects, we must distinguish the two type of one-loop corrections, 
corresponding to I\{x) and I\{x + £) respectively; in the latter case the factor is in the integrand 
of the / integration. We now demonstrate that neither type of contribution contains collinear 
divergences. 

We start by considering the terms proportional to I\{x) arising from diagrams with a gluon 
attached to the outgoing quark line with momentum xq^. The corresponding diagrams are (12), 
(2V)' and (27). Repeating the steps in section 4.1.1 we find 

2 h{x) f dH x-f 1 



(12) - (2V)' = -VC, J — _> _ ^—^ . (58 ) 

The remaining contribution of this type comes from diagrams (27). The numerator of the light- 
quark propagator has a factor of 2{x — C)Q2 — ^(^ — ^> wriere is the Lorentz index of the 
gluon with momentum I. Consider first the diagram in fig. 8 (a), for which I = k% — k±. Using 
the Ward identity 
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inside the quark loop, we readily see that the integral over the quark loop can be reduced to 
ones with only three quark propagators. The first term on the right-hand side of eq. (59) gives 
a contribution proportional to I± (x) . Similarly consider the insertion of the gluon as in fig. 8 (b) 
and the Ward identity 

T^— ~ h) T^— = T~ T~~ — • ( 6 °) 

f 3 — m f 4 — m f 3 — m f 4 — m 

The first term on the right-hand side of eq. (60) cancels the second term of eq. (59), and the 
second term of eq. (60) gives a contribution proportional to I\(x + £). These arguments can be 
generalized to include all the insertions of the photon and gluon onto the quark loop. The term 
proportional to h(x) precisely cancels the contribution from (12) — (2V)' in eq. (58). 

The cancellation of the terms with integrands proportional to I\ (x + £) in diagrams with the 
gluon with momentum I attached to the external light quark with momentum xq 2 proceeds in 
a similar way. The corresponding contributions are as follows: 

(23) -(2V)' = ig*C A r dH *-^ + ^ 1 



(26) = ig 2 C A 



(2vr) 4 £ x + £ I 2 (xq 2 - I) 2 ' 
dH x-C h{x + 1 1 



(2vr) 4 x x + £ I 2 (xq 2 - I) 2 ' 
(07\ in 2 r f dH x-g /xQr + Q l 1 

(27) = - l9CA )^f- X I* T7)2 



so that 



(23) - (2V)' + (26) + (27) = (61) 
and all the mass singularities cancel as required. 

The cancellation of the mass singularities from the diagrams in which the gluon with momentum 
I is attached to the light antiquark with momentum xq 2 proceeds in a very similar way, except 
that now we have factors I\{x) and I\{x — £). For the terms with a factor I\{x) the contributions 
are: 

2ri h{x) f dH x-i 1 1 



(34 ) - (4F)< = -#C A ^-1J—^ T2 T — f? , (62) 



... h(x) f dH 1 1 



(46) = - u/ c A -^J— pJ — J? , (63) 
fA*\ ■ 2/~< h(x) f dH x 1 1 

,47) = ^_Uy^ ?? ^_ 7p , (64) 

(65) 

which sum to zero, and for the terms with a factor I\(x — £) in the integrand we have 

< M > " (4F) ' " - (47) " I §? ^ ? T^W ' (66) 

Thus again, all the collinear divergences from the region in which a gluon's momentum is parallel 
to q 2 cancel. 
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4.1.3 Divergences from the Region Collinear to k 



For the kinematical situation considered in ref. [6], in which we keep m^/m 2 , but neglect 
<72 • k/m 2 , the cancellation of the mass singularities from the region in which a gluon is collinear 
to k follows exactly as in the case of the operator Qg discussed in section 4.1.3 (see equations 
(45) -(48)). 

This concludes the demonstration of the cancellation of mass singularities. However it is also 
instructive to ask whether they would have cancelled if we had not neglected qi ■ k/ m 2 . We now 
rewrite the amplitude at lowest order as 

A? = ^-I 2 (x,k + ), (67) 

where I%(x, k + ) = k + Ii(x, k + ) and the k + dependence in eq. (57) and subsequent equations was 
implicit. For all the discussion above, the k + dependence in h(x) was contained in a simple 
overall factor of l/k+. Now, if we keep terms of 0{q 2 • k/m?), we have two types of term in 
the integrands of the loop integrations, those proportional to l2(x,k+) and those proportional 
to I 2 (x, (1 — cr)k + ). We now demonstrate that the corresponding mass singularities cancel also 
in this case. 

The distribution amplitude of the initial state is defined by 

*%{k+) = J dz- e fk + z - (0\<? p (z)[z,0}b a (0)\b(p)<?(k))\ z+ , z±=0 , (68) 
so that at tree level (see eq. (18)) 

<S>% { °\k + ) = ^5{k + - k + )vp{k)u a {p-k) , (69) 

where v and u are the free spinor wave functions of the light antiquark and the 6-quark respec- 
tively. In the collinear region which we are considering here, the one-loop correction from the 
diagram in fig. 5(c) can be written in the form 

Sggfo.) = 2ig 2 C F v f3 (k)u a (p) J -^L ii-^ { 6(k + - k+) - 5(k + - (1 - a)k+) } , 

(70) 

where the subscript (c) denotes that this is the contribution from the diagram of fig. 5(c), and 
the superscript (1) that it is a one-loop contribution. W <g) <g) (°) is therefore 

(3B)' = 2ig 2 C F -L J \ { (1 - a)I 2 (x, k+) - I 2 (x, (1 - a)k + ) } • (71) 

The results in section 3.2.2 are now modified as follows: 

(23) - -2 % g [C F -— J—/^ ~ 72 {k _ l? > ( 72 ) 





1 1 




xk + 




1 1 


2 J 


xk+ 



(34)= 2^[ Cf -^|^//LI^|£W (73) 



... i r an h{x,k+) 



< 36 > = -« c «^Jwf(W=W- < 74 > 
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Finally we have to consider the diagrams (37) which no longer vanish when we include the terms 
of 0(q2 • k/ml). Using the Ward identity in eq. (59) the contribution from these diagrams is 
readily found to be: 

. 2 r 1 f dH I 2 (x,k + )-I 2 (x,(l-a)k + ) 1 

Thus we have 

(23) + (34) + (36) + (37) = (W)' , (76) 
so that there are no mass singularities remaining in the hard-scattering amplitude. 



5 Conclusions 

In this paper, we have studied the radiative decays B — > Vj (where V = p, K*) in the framework 
of QCD factorisation. We focused on spectator interactions, where the most significant contri- 
butions come from the chromomagnetic operator Qs and the four-quark operator Q±. These 
contributions had previously been considered explicitly at leading order [6, 7], and they exhibit 
similar features to the purely radiative decays B — > j£v (77, r y£ + £~~) that we have already con- 
sidered [11, 12]. However, the presence of two light-cone distribution amplitudes (those of the 
B and vector mesons) leads to technical and conceptual modifications. 

We have performed an explicit next-to-leading-order computation of the contribution from Q$ 
in the heavy-quark limit mi, — ► 00 (specifically we have calculated all the terms containing mass 
singularities and large logarithms), and we showed that the spectator interactions factorise, i.e. 
they can be written as the convolution of a hard-scattering kernel, computable in perturbation 
theory, and of two light-cone distribution amplitudes (one for each meson) describing soft physics. 
The explicit results are presented in appendix A. By studying the soft and collinear regions of 
the loop momentum, we can understand the factorisation of mass singularities using a heuristic 
argument (presented in section 3). The same argument can be applied to the 4-quark operator 
Qi, so that its contribution to spectator interactions are also expected to factorise at next-to- 
leading order. In appendix B we use the collinear Ward identity to demonstrate the factorisation 
of mass singularities at next-to-leading order for all the weak operators contributing to B — > Vj 
decays. 

This success of QCD factorisation for B — > spectator interactions leaves several questions 
open. First, the presence of an intermediate scale a/ ' A.QCD m bi corresponding to the virtuality 
of the exchanged gluon, yields large (Sudakov) logarithms at all orders of perturbation theory. 
The latter could be resummed using renormalization-group arguments within the Soft-Collinear 
Effective Theory [13, 14, 8]; for the resummation in purely radiative decays see refs. [11, 12]. 
However, the mixing of the operators and hence the subsequent resummation are much more 
intricate than in the case of b — ► u transitions 5 . We have only considered twist- 2 distribution 
amplitudes for both mesons: higher-twist distribution amplitudes deserve a separate study. 
One could also consider the interesting decay B — > V£ + £~, where the light vector meson is 
longitudinally polarized and the lepton pair has a small invariant mass. In ref. [7], the spectator 
interactions were analysed at leading order: they receive contributions from additional diagrams 
to those considered here, involving (in particular) the second leading-twist i?-meson distribution 
amplitude Of course, one would also like to understand in detail hard spectator interactions 
in two-body nonleptonic -B-decays. 

5 We thank M. Neubert for instructive discussions on this point. 
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A Contribution of Qg at one loop 

In this appendix we present the results from explicit evaluation of the Feynman diagrams at one- 
loop order for the chromomagnetic operator Q$. Specifically we exhibit the mass singularities 
and all the large logarithms. We denote as (ij) the contribution from the diagram where i 
and j are connected with a gluon in Fig. 3. Each diagram's contribution is given in units of 
a s /(4vr)^ 0) . We use dimension regularisation to regulate both the ultraviolet and infrared 
divergences, but introduce separate scales for the two cases. 

A.l Abelian Component 

We start by considering those diagrams which have an Abelian component, i.e. which have a 
component proportional to Cf- 



A. 1.1 Contributions to 

There are four diagrams with a gluon which links an initial-state quark with a final-state one: 

C A 



(12) 

(14) 

(23) 
(34) 



C F 



C F - — 



C f - 



C a 



Cf- — 



Mir 



1, o Ml , , MV 
--log 2 ^-21ogxlog^ 

2 Mir 

1 , o Ml 

2 /'ft< 



2 , 2(k -q 2 ) 2 2 ( k ■ Q2) 
o lo S -^ + (21ogx-2)log— -2 log m2 

9 2(/c ■ q 2 ) 
log 2 \ 2 7 + 2 log x log 



Mir 
2{k ■ q 2 



Mir 
log 2 2{k • 92 



Mir 

01 ^A\^ 2(fc • q 2 ) , 2(fc • g2 ) 
+ (-2 log x + 4) log § log g 



Mir 



Mir 



Mtjv 



(77) 
(78) 
(79) 
(80) 



One may also link the external lines of the same hadron, yielding 

2\ 



(13) = C F 
(24) = 



-iiog 2 -f^y 

2 S \ \M Bf i 1R J 



at leading twist. 
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Figure 9: An additional diagram, D, which contains a large logarithm and which is not included 
in the set of diagrams { (ij) } defined in section 3. 



When an external line is linked to the internal quark propagator one gets 



(15) 
(25) 
(35) 
(45) 





C F 


C F \og 



, M B 

21og^ 
/% 

2(fc • q 2 ) 
Mir 



i M l 
Muv 



2ilogx + 2 

x 



There are also contributions from the wave function renormalization of the external lines 

(11) = C F 



1, Ml , Ml 
2 lo gT^ + lo g^ 



(22) 
(33) 
(44) 
(55) 



= 



Mtjv 
C F 



Mir 



C f 
C f 



1 Mtjv 

2 Muv 

1 1q Mir 

2 Muv 



Muv 



where the vanishing values correspond to massless tadpole diagrams in dimensional regularisa- 
tion. 

Finally, there is a contribution containing large logarithms that cannot be obtained by adding 
one gluon to the leading- twist tree- level diagram and is therefore not included in the set { (ij) } 
defined in section 3. We denote it by D and exhibit it in fig. 9; its contribution to the amplitude 
is: 



D = (C F - C a /2)- 



-2 log 



2k -q 2 

Ml 



log 



M 



Muv/ 



(81) 



The other diagrams not included in the set { (ij) } do not give large logarithms when ^uv is set 
equal to Mg. 



A. 1.2 Distribution amplitudes 

We consider now the contribution to the distribution amplitudes convoluted with the lowest- 
order kernel, W <g> T 8 {0) ® (°) + (°) ® T 8 (0) <g> <W (l) . As explained in section 3, we 
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denote the diagrams in a similar way to those for the amplitude. Once again we use dimensional 
regularisation to regulate both infared and ultraviolet divergences. The infrared scale, fiiR, is the 
same for the matrix elements and the distribution amplitudes. The ultraviolet scale is different 
in the two cases however: for the matrix elements ^uv corresponds to the renormalization scale 
(Ir, whereas for the distribution amplitudes z^uv corresponds to the factorisation scale [i F (see 
ref. [11] for further discussion). In agreement with our definition of & H , we use the HQET 
Lagrangian when a gluon is connected to the 6-quark. 

We can link an external line to the path-ordered exponential of one distribution amplitude 



(IB)' 

(2vy 

(Avy 



c F 



1 2 2fc 2 i 
"2 l0g ~ 



= C F 
= C F 



"uv 
2 log 

Mir. 
Mir 



2 log 



uv 



C F log — r 
Mir 



We can also link two external lines 

(13)' = C F 
(24)' = 



--log 2 

2 & 



2-logx + 2 

x 



f 2k-q 2 V 



There are also contributions due to the wave function renormalization for external lines 



Mir 



(11)' = C F 

(22)' = = C F 

(33)' = = C F 

(44)' = = C F 



z ^uvJ 
Mm 



log 



v 



uvJ 

log^' 
^UVJ 



where the vanishing values correspond to massless tadpole diagrams in dimensional regularisa- 
tion. 

We could also try to link two points of the same path-ordered exponential. These diagrams 
involve the propagator associated with A + (az)A + ((3z) or A~(az)A~ ((3z), which vanishes in the 
Feynman gauge. 



A. 2 Remaining Nonabelian Components 

We now evaluate the one-loop contributions to „4g^ whose colour factor is proportional to Ca- 
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A. 2.1 Diagrams with one Three-Gluon Vertex 



The contributions from diagrams in which a gluon links an external quark line to the gluon one 
through a three-gluon vertex are as follows: 
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-2 l0g " 
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M _, + log , jlog ^ 
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3. 2(fc- g2 ) , 01 2(fc • g2 ) 
~o !og— 2 +2 log g 



2 Mfjv ^IR 
keeping in mind that (46) is the same diagram as (36). 
In the case of the internal quark propagator, we get 

(56) = C A 



1 

2x 



3 l + x 

4 ~ 



Muv- 



n\ Mr 

log^ 
Muv- 



A.2.2 Diagrams for the Gluonic Vacuum Polarization 

The gluonic wave function renormalization yields 

2{k ■ q 2 ) 



(66) 



log- 



Muv 



A. 2. 3 Diagrams with two Gluons attached to Qg 

We can also attach two gluons to the vertex with the operator Q§. One of them interacts with 
the spectator quark, whereas the second is attached elsewhere. 

When one of the gluons is attached to an external line, we obtain 



(17) 
(27) 
(37) 
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where (37) and (47) denote the same diagram. 

The result for diagram (57) obtained by attaching one of the gluons to the internal quark 
propagator is 



(57) =C A x— log 



Mr 



4x Muv 

and for diagram (67) in which the gluon is attached the internal gluon propagator is 

^ 3 , 2(k ■ q 2 ) 
(67) = C A x - log 2 . 

z Muv 
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A. 3 One-loop results 



We now combine the results presented above to obtain the large logarithms at NLO in the 
amplitude, in <j? B <g) T§ <8> 3>j_ and hence finally in the hard-scattering kernel. These next-to- 
leading-order results contain logarithms of ratios of the scales (j,r, fj,p, Mr and q2-k. Following 
ref. [11, 12], the natural choices for the renormalisation and factorisation scales are 



jji R = 0(M B ) Hf = m = 0(q 2 ■ k) 



(82) 



But even with this choice, large logarithms will remain due to the presence of three distinct scales. 
Therefore, in presenting our results we focus on the mass singularities and large logarithms. The 
amplitude up to one- loop order is: 
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where (//r) denotes the tree- level matrix element with /j,r as the scale of the strong coupling 
constant, and the ellipses denote terms without large logarithms. 

The contribution from the distribution amplitudes is 
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(83) 
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The one-loop result for the hard-scattering kernel at next-to-leading order is therefore: 

.(o+i) „ ^ _ T ^\k + ,u- m = ^)x 

5 



(84) 
+ 0(a 2 ). 
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where Tg°\k + ,u; [1r) denotes the tree-level hard-scattering kernel with [1r as the scale of the 
strong coupling constant, and the ellipses denote terms without large logarithms. The key point 
is, of course, that there are no mass singularities in the expression for T$. 



B Cancellation of Collinear Divergences and the Collinear Ward 
Identity 

In this appendix we show that the cancellations of collinear divergences demonstrated in sec- 
tions 3 and 4 are a general consequence of the collinear Ward identity and hold for all the 
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Figure 10: Amplitude for the auxiliary process bq — > qq^f + gluon.. The spring represents the 
photon (with momentum q\) and the curly line the gluon (with outgoing momentum /). The 
first four diagrams are those in which the final-state gluon is attached to one of the four external 
quark lines. The fifth diagram represents all the remaining ways in which the final-state gluon 
can be attached. 



weak operators contributing to B — ► V7 decays 6 . It will be useful to consider (in perturbation 
theory) the on-shell amplitude bq' — > qq ,r y + gluon, where the momentum of the gluon is /. We 
represent this amplitude by the graphs in Fig. 10, exhibiting explicitly the diagrams in which 
the gluon is attached to an external quark line. We write this amplitude as £g(l)Ai^, where e g 
is the polarization vector of the gluon. The Ward Identity we wish to exploit is 

l»M f , = 0. (86) 

The process we are actually studying in the evaluation of the one-loop contributions to the hard- 
scattering kernels is bq' — > qq'j, where the momenta of the incoming 6-quark and light antiquark 
are q\ + q2 — k and k respectively and those of the outgoing quark and antiquark are xq2 and 
xq2- We denote this amplitude at lowest order in perturbation theory by A^°'(x, k), without 
specifying which operator mediates the weak transition. The corresponding hard-scattering 
kernel is written as T®\ where we suppress the spinor indices. 

We start by the demonstration of the cancellation of collinear divergences from the region in 
which / is collinear to q2, I — iq2- It is convenient to use the identity in eq. (59) to represent 
the Ward identity (86) by Fig. 11. On the right-hand side we have eliminated the fermion 
propagator immediately to the left of the solid dot. Since the external particles are on-shell we 
can use the identity in Fig. 11 with the external wave functions amputated. 

We now consider the identity in Fig. 11 with the external wave functions amputated, and with 
the momentum of the final state quark equal to xq2 — I. We multiply both sides of Fig. 11 by 
the factor 

-*°^r vh^w ' (87) 

6 Similar arguments were presented in a different context in Ref. [23]. 
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— I X 




Figure 11: Representation of the Ward identity in eq. (86). The solid dot in the two diagrams 
on the right-hand side indicates that the gluon is attached with factor igT a (but no 7-matrix) 
and the propagator to the left of the dot has been eliminated using eq. (59). 



and the outgoing quark line by the colour matrix T a (a is the colour label of the final-state 
gluon). Integrating over I (in the collinear region), on the left hand side we generate all the 
diagrams with a gluon attached to final-state quark line (i.e. the line denoted by 2), except 
for (22) and (24) 7 . From the identity in Fig. 11, we then immediately understand that for 
any operator the sum of these diagrams gives a contribution which can be written in terms of 
a convolution of the properties of the final state with the lowest-order amplitude. Specifically, 
Fig 11 gives 



n^2,4 



2ig 2 C F 



d 4 l 
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(2vy, 



(89) 
(90) 



where the second line follows from eq. (36). Here we have restored the external quark wave 
functions. 



Following similar arguments and eq. (37) one also obtains 

dH x-£ A(°\x,k) -A(°\x 



£(4n) 



2ig z C F 



(2tt) 4 £ 



l 2 (xq 2 - If 
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(91) 

(92) 
(93) 



Finally one can use the same procedure for the collinear divergences from the region in which I 
is parallel to k (I ~ ak). In this case we use the identity in eq. (59) to eliminate the propagators 

7 Apart from the factor in eq. (87), we have set (xq2 — I) 2 = everywhere, which is correct at leading twist. 
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to the right of the gluon vertex in the first two diagrams of Fig. 10, and obtain 



(3n) = 2ig 2 C F 

n^l,3 




(94) 



$g (1) ®T(S$ w "' (0) 
(3B)> . 



(95) 
(96) 



Thus all the collinear divergences are absorbed into the light-cone distribution amplitudes for 
any of the weak operators Qi (the argument presented here considers all singularities from 
the collinear region, and includes in particular the soft-collinear modes discussed recently in 
refs. [24]). It is straightforward to verify that for Q 8 , for which 



one recovers the results in section 3 and similarly for Q\ one recovers the results in 4. 
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